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Context. An accelerated logarithmic potential models the mean motion of stars in a flat rotation curve galaxy that 
sustains a wind system. For stars outside the galactic wind launching region, the asymmetric removal of linear momen- 
tum by the wind is seen as a perturbing acceleration superimposed onto the galactic potential. 
Aims. We study the integrability of stellar motion in an accelerated logarithmic potential. 
Methods. We use surfaces of section of the dynamical system to probe the integrability of motion. 
Results. We provide numerical evidence that motion in an accelerated logarithmic potential is non-integrable. 
Conclusions. Large scale chaotic diffusion occurs for lower values of the projected angular momentum along the direction 
of acceleration and persists at all values of the angular momentum in the outer part of the galaxy inside the truncation 
radius where the galactic acceleration balances the wind-induced acceleration. 
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1. Introduction 

The advent of new observations of galactic winds have 
shown that these processes are ubiquitous features in star- 
burst and active galaxies (Veillet et al. 2005). In this paper, 
we lay down the dynamical premise for studying the ef- 
fect of asymmetric linear momentum removal from a galaxy 
with a flat rotation curve. Our investigation is motivated 
by the recent studies of the dynamical influence of asym- 
metric stellar jets on protoplanetary disks and the forming 
planets they contain (Namouni 2005, 2007; Namouni and 
Zhou 2006; Namouni and Guzzo 2007). These works have 
shown that the jet-induced acceleration leads to a number 
of dynamical consequences such as strongly modifying pro- 
toplanetary disk structure profiles, forcing the radial con- 
traction and expansion of the disk, and enhancing the inter- 
nal heating of the disk as it follows its varying state of least 
energy induced by the jets' variablity. Jet-induced acceler- 
ation also helps truncate the protoplanetary disk, excites 
the eccentricity of forming planets, and induces radial mi- 
gration in binary systems. In view of the variety of these 
consequences, we expect that the wind-induced accelera- 
tion alter star orbits and affect the collective response of 
the galactic disk. 

Modeling the motion of a star in a flat rotation curve 
galaxy as that of a test particle moving in a central loga- 
rithmic potential, logr, where r is the distance to the galac- 
tic center (Binney and Tremaine 1987), the asymmetric re- 
moval of linear momentum from the galaxy can in turn be 
modeled by an additional acceleration that does not depend 
on the position and velocity of the star. This approximation 
is valid as long as we consider the motion of star outside 
the wind launching region. In effect, as such a star revolves 
around the matter inside its orbit, it is only affected by 



the variation of mass and momentum losses inside its or- 
bit. It is worth noting that the effect under consideration 
is not that of the direct action of the galactic wind on the 
star such as that of the wind blowing on the star -which is 
physically negligible. Rather it is the indirect effect of the 
loss of linear momentum by the galactic disk on the motion 
of stars outside the wind launching region. 

To set up the dynamical premise of this process, we first 
discuss the zero- velocity curves of the accelerated logarith- 
mic potential as well as its orbits of least energy (section 2). 
The surfaces of section allow us to show that orbital motion 
in the accelerated logarithmic potential is not integrable 
(section 3) in contrast to the accelerated Kepler poten- 
tial (Epstein 1916, Sommerfeld 1929, Landau and Lifschitz 
1969) that models orbital motion in the presence of a stellar 
jet (Namouni and Guzzo 2007) . Section 4 contains conclud- 
ing remarks. 

2. Zero velocity curves and least energy orbits 

Galaxies with asymmetric wind systems lose linear momen- 
tum in the region inside the wind launching radius. The net 
momentum loss accelerates the wind launching region with 
respect to the matter outside it in the direction opposite to 
momentum loss (Namouni 2005, 2007; Namouni and Guzzo 
2007) . In essence, asymmetric momentum loss from a wind 
system accelerates the inner part of the galaxy like a rocket 
with respect to its outer part. Modeling the gravity of a flat 
rotation curve galaxy as a logarithmic potential (Binney 
and Tremaine 1987), the motion of a star outside the wind 
launching region follows the equations: 



d x v~ 
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Fig. 1. Zero velocity curves for the value of h 2 z 
0.058313294. 
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Fig. 2. Sphere of least energy orbits. The solid (dashed) 
line denotes stable (unstable) orbits. 



where x denotes the position of the star in the galaxy, v c 
is the galactic circular velocity and A is the wind-induced 
acceleration. We examine the effect of momentum loss dur- 
ing a stationary wind episode and therefore take A to be 
constant in time. The presence of the acceleration A intro- 
duces natural units for time and distance. These are given 
as: 



V ~T ^d T=% 
A A 



(2) 



The former is related to the boundary of the logarithmic 
potential set by the equality of the gravitational accelera- 
tion and the wind-induced acceleration. The latter is the 
typical excitation time for stellar orbits in the presence of 
acceleration. In terms of these units, the equations of mo- 
tion are given as: 



ci"x 
ch2~ 



(3) 



where e z is the unit vector along the direction of accelera- 
tion A that we choose as the z-axis of the inertial reference 
frame. Note that it is still possible to keep track of the rela- 
tive strengths of the logarithmic potential and the constant 



acceleration. For instance the logarithmic potential is dom- 
inant for |x| <C 1, and a unit orbital period corresponds to 
motion near the boundary where the gravitational acceler- 
ation balances A. 

As the acceleration is constant, the corresponding force 
derives from the potential R = z and hence the dynamical 
system ([3]) admits the constant of motion: 



p 2 9 2 + z 2 



\ log (P 2 



(4) 



where p, 9 and z are the usual cylindrical coordinates. In 
addition to the energy E, the projection of the angular 
momentum h = x x v along the direction of acceleration 
(or equivalently the z-axis), h z — h e z = p 2 6, is a constant 
of motion. Substituting the expression of 9 in terms of h z 
into Equation (T?|), the energy equation now reads: 
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The zero velocity curves of the system ([3]) are obtained 
from the effective potential of the energy equation - the 
last three terms of (JS|). In Figure (1), we show such curves 
for the value of h 2 z = 0.058313294. The region of bounded 
motion is restricted in space inside a radius of p 2 + z 2 ~ 1 
where the gravitational attraction exceeds the amplitude 
of the perturbing acceleration. There are two equilibrium 
(least energy) points located at (p = 0.25, z = 0.067) and 
(p = 0.452, z — 0.714). The former is a stable orbit and 
the latter is unstable. 

We determine the properties of least energy orbits by 
equating with zero the partial derivatives of E with respect 
to p, z,p and z. The derivatives with respect to p and i 
show that least energy orbits are circular and planar. The 
derivative with respect to p gives the expression of h z : 



hi = 



p 2 + z 2 



(6) 



or equivalently the orbital rotation rate fl = 9 = h z p~ 2 in 
terms of p and z. Lastly, the derivative with respect to z, 
which is the vertical projection of the equation of motion, 
shows that least energy orbits hover above the center of 
gravitational attraction according to the equations: 



l±y/T^4p 2 



(7) 



From this equation, it is seen that the locus of least energy 
circular orbits is a sphere centered on (p = 0, z = 1/2) 
of radius 1/2 (Figure 2). The upper and lower signs do 
not correspond to the stability of the equilibrium orbits. 
To determine which orbits are stable on the least energy 
sphere, we first note that every least energy orbit shown 
in Figure (2) has a different value of the projected angular 
momentum, h z , and the energy, E. From the zero velocity 
curves (Figure 1), we know that for a given value of h z there 
are two least energy orbits, one stable and one unstable — 
the smallest value of h z , zero, must corresponds to the two 
least energy points at p — and z = (stable) and z = 1 
(unstable). Using Equation (J7J), it is possible to express h z 
Q as a function of z for least energy orbits and solve for 
the values of z (and hence p) that determine the stable and 
unstable orbits. This substitution yields: 



h 2 z = z{z-lf 



(8) 
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Fig. 3. Surfaces of section for h\ = 0.058313294. The energy levels are E = -0.75 (top left panel), E = -0.80 (top right), 
E = —0.90 (bottom left) and E = —0.95 (bottom right). The solid line is the sombrero profile z = (1 — \J\ — 4p 2 )/2 
that divides space into z < orbits (above the line) and z > orbits (below the line). 



This shows that h 2 z has a maximum at: 

z h = 1/3, and p b = v/2/3, (9) 

and as expected h z vanishes for z = and z = 1 (fi = 0). 
The maximum value of ft. 2 is 4/27 and corresponds to the 
boundary of stable least energy orbits. It also defines the 
truncation radius of the galaxy. The condition 

|/i 2 |<2/3v / 3 (10) 

is therefore necessary for bounded motion to occur. The 
profile of stable and unstable orbits is shown in Figure (2) 
with solid and dashed lines respectively. The locus of the 
stable orbits in the pz-plane has a profile similar to (but 
less flat at its base than) the sombrero-shaped profile asso- 
ciated with the accelerated Keplerian potential (Namouni 
2007; Namouni and Guzzo 2007). For consistency with the 
Keplerian case, we call the locus of the stable orbits of an 
accelerated galactic potential, the sombrero profile. An ad- 
ditional difference between the accelerated Keplerian and 
logarithmic potentials is that the former admits a global 
third integral of motion implying that orbits are everywhere 
regular -even in the vicinity of the zero-velocity separatrix 
and the unstable least energy orbit. The accelerated loga- 
rithmic potential does not have a global third integral. This 
we show in the next section. 

3. Surfaces of Section 

The integrability of stellar motion in the accelerated loga- 
rithmic potential may be determined numerically by draw- 



ing surfaces of section of the dynamical system . Sections 
in phase space are chosen for (i) a constant energy, E, (ii) 
a constant projected angular momentum, h z , and (iii) with 
the condition z — 0. Sections usually require a crossing di- 
rection such as that given by z > 0. However, as we find 
it convenient to represent sections in the pz-plane, we can 
readily determine the sign of z everywhere in that plane. As 
z equals the vertical component of the gravitational poten- 
tial augmented by the perturbing acceleration, z vanishes 
exactly on the sphere of least energy (Equation[7]). The sign 
of z is negative inside the least energy sphere and positive 
outside it. This allows us to explore numerically the whole 
area inside the zero velocity curve associated with the val- 
ues of E and h z , and know the sign of z. Accordingly, we 
plot the curve of z = (Equation [7]) on all surfaces of sec- 
tion which are now allowed to include both z > orbits 
and z < orbits. 

Using the same value of the sombrero orbit as that of 
Figure (1), p = 0.25, we set h\ = 0.058313294 and compute 
the surfaces of section for E = —0.95, —0.90, —0.80 and 
—0.75. The value E = —0.95 corresponds to a small region 
around the sombrero orbit. Increasing the value of E widens 
the area of possible motion. Orbits have been computed for 
a duration t = 200 as well as for various values of p, z and 
z whereas p is obtained from the energy equation. 

The section for E = —0.75 (Figure 3, top left) shows 
that motion in the accelerated logarithmic potential is not 
integrable with large chaotic regions developing for smaller 
and larger elevations z. Irregular motion is milder near the 
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Fig. 4. Surfaces of section for different values h 2 — 6.24 x 10 -4 (top row) and h 2 = 0.128 (bottom row). For h 2 = 
6.24 x 1(T 4 , the energy levels are E = —0.97485 (top left) and E = —1.5 (top right). For h 2 z = 0.128, the energy levels are 
E = -0.590575 (bottom left) and E = -0.596 (bottom right). The solid line is the sombrero profile z = (1- y/l - V)/2. 



plane z = that contains the gravitational center. Noting 
that the surface of section is asymmetric with respect to 
the sombrero profile z = (1 — ^/l — 4p 2 )/2, restricting the 
surface of section to positive z only would have lead to a 
loss of information on the motion inside the least energy 
sphere. The lowering of the value of E shrinks down the 
region of possible motion and as well as that of irregular 
motion (Figure 3, top right). In particular for E = —0.90 
motion is further restricted to a small region around the 
sombrero orbit (p = 0.25, z = 0.067) and an unstable point 
appears very near it (Figure 3, bottom panels). Lowering 
the energy further makes the unstable point converge to- 
wards the sombrero orbit (E = —0.95, Figure 3, bottom 
right) implying the presence of a mild separatrix-type dif- 
fusion in the vicinity of the least energy orbit. 

The extent of the regions of chaotic motion also depends 
on the value of the angular momentum h z . We find that 
the smaller the projected angular momentum h z the larger 
the chaotic motion regions. This is shown in Figure (4) 
where we plot surfaces of section for h 2 = 6.24 x 10~ 4 
and 0.128. The former corresponds to a stable sombrero 
orbit at p = 0.025 and z = 6.25 x 10~ 4 and the latter 
to p = 0.4 and z — 0.2. On Figure (1), these two orbits 
are on each side of that corresponding to Figure (3). For 
h 2 z = 6.24 x 10~ 4 , the energy levels are E = -0.97485 (top 
left) and E = —1.5 (top right). The former has been chosen 



close to the energy value of the unstable least energy orbit. 
The region of chaotic motion is large and remains so even 
when the energy value is lowered (E = —1.5) so that the 
zero velocity curve extent (i.e. the size of the surface of 
section) is 4 times smaller in the vertical direction that it 
is for E = -0.97485. For h 2 z = 0.128, the energy levels 
are E = -0.590575 (bottom left) and E = -0.596 (bottom 
right). As for the previous value of h z , the former energy has 
been chosen close to the value of the unstable least energy 
orbit. This time, the region of chaotic motion is confined to 
the vicinity of the unstable orbit and disappears when the 
energy value is lowered slightly (E — —0.596) so that the 
zero velocity curve extent is only 1.4 times smaller in the 
vertical direction that it is for E = —0.590575. We conclude 
that for the smaller h z , the extent of the chaotic region is 
much larger than that of Figure (3) and irregular motion 
dominates the surface of section even when the energy value 
is far from that of the unstable least energy orbit. For the 
larger value of h z , extended chaotic motion is confined to 
the vicinity of the unstable least energy orbit as soon as the 
energy level lowered away from that of the unstable orbit, 
no region inside the zero velocity curve shows large chaotic 
motion. 
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4. Concluding remarks 

In this paper we showed that stellar motion in an acceler- 
ated logarithmic potential is not integrable unlike motion 
in an accelerated Keplerian potential. We have used zero- 
velocity curves, least energy orbits as well as surfaces of 
section to determine the global topology of orbital motion. 
We find that for a given projected angular momentum, the 
larger the area of the zero-velocity curve in pz— space, the 
larger the fractional size of the chaotic motion region. The 
size of the chaotic motion regions also depends on the value 
of the projected angular momentum h z . The smaller h z , the 
larger the extent the chaotic motion region. Smaller values 
of h z correspond to inner sombrero orbit (i.e. p -C 1) while 
larger values of h z correspond to zero velocity curves whose 
corresponding motion is confined near to the truncation ra- 
dius of the galaxy (i.e. the stable and unstable least energy 
orbits become close). 

For the diffusion that we have uncovered in this work 
to affect stellar orbits, the local diffusion time has to be 
smaller than the age of the galaxy. This implies that in 
the outer regions where the dynamical time (i.e. the orbital 
period 2itr/v c ) is comparable to the age of the galaxy, dif- 
fusion does not have enough time to occur. The reason is 
that the local diffusion time is always and everywhere larger 
than the local dynamical time. For instance, the scatter of 
points in the top left panel Figure (3) has been obtained 
after 200 time units -a time unit being the dynamical time 
at the truncation radius of the galaxy, T ©. For a typi- 
cal galaxy, such a timespan is larger than the age of the 
Universe. 

Diffusion therefore affects stellar orbits in an accelerated 
flat rotation curve galaxy for low values the projected angu- 
lar momentum h z and well inside the truncation radius (i.e. 
the outer edge of the sombrero profile) . Stellar orbits near 
sombrero orbits are more regular as diffusion is negligible. 
In these parts of the galaxy, it is possible to make use of 
secular perturbation techniques to ascertain the evolution 
of stellar orbits whose periods are small compared to the 
duration of the wind episode. 

Acknowledgements. The authors thank the referee for useful com- 
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